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ABSTRACT 

The iPTF detection of the most recent outburst of the recurrent nova system RX J0045. 4+4154 in the An- 
dromeda Galaxy has enabled the unprecedented study of a massive ( M >1.3 M 0 ) accreting white dwarf (WD). 

We detected this nova as part of the near daily iPTF monitoring of M3 1 to a depth of R « 21 mag and triggered 
optical photometry, spectroscopy and soft X-ray monitoring of the outburst. Peaking at an absolute magnitude 
of Mr = -6.6 mag, and with a decay time of 1 mag per day, it is a faint and very fast nova. It shows optical emis- 
sion lines of He/N and expansion velocities of 1900 to 2600 km s _1 1-4 days after the optical peak. The Swift 
monitoring of the X-ray evolution revealed a supersoft source (SSS) with kT e g « 90- 110 eV that appeared 
within 5 days after the optical peak, and lasted only 12 days. Most remarkably, this is not the first event from 
this system, rather it is a recurrent nova with a time between outbursts of approximately 1 year, the shortest 
known. Recurrent X-ray emission from this binary was detected by ROSAT in 1992 and 1993, and the source 
was well characterized as a M > 1.3 M 0 WD SSS. Based on the observed recurrence time between different 
outbursts, the duration and effective temperature of the SS phase, MESA models of accreting WDs allow us to 
constrain the accretion rate to M > 1 .7 x 1CT 7 M 0 yr _1 and WD mass > 1 .30 M 0 . If the WD keeps 30% of the 
accreted material, it will take less than a Myr to reach core densities high enough for carbon ignition (if made 
of C/O) or electron capture (if made of O/Ne) to end the binary evolution. 

Subject headings: galaxies: individual (M31) — novae, cataclysmic variables — supernovae: general — X- 
rays: binaries — white dwarfs 


1. INTRODUCTION 

Classical novae are the observable outcome of unstable 
thermonuclear burning on accreting white dwarfs (WDs). 
In most tight binaries with WDs of typical masses of M = 
0.8 Mq, the recurrence time is tens of thousands of years. 
However, if the accretion rate, M, is high and the WD mass is 
large, the time between flashes can become short enough so 
that the recurrence can be measured. Due to the large accre- 
tion rate and insignificant ejection mass loss, recurrent novae 
(RNe) have been proposed to grow toward the Chandrasekhar 
limit and could be promising progenitors of Type la super- 
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novae (Starrfield et al. 1988; Di Stefano 2010). Within our 
galaxy, these recurrent novae (RNe) have, at the minimum, a 
recurrence time of 10 years (Schaefer 2010). 

As the nearest large galaxy neighbor of the Milky Way, 
the Andromeda galaxy (M31) provides the best opportunity 
for studies of classical and recurrent novae. Extensive pho- 
tometric and spectroscopic surveys have been conducted to 
search for novae in M31, and have resulted in the discovery 
of over 900 novae, with over 100 having well-sampled light 
curves, optical spectra, or X-ray observations (Hubble 1929; 
Arp 1956; Darnley et al. 2004; Henze et al. 2010; Shafter 
et al. 2011; Cao et al. 2012, see also the website maintained 
by Pietsch 14 and references therein). There are six confirmed 
RNe in M31, and a few other strong candidates (Shafter et al. 
2013). 

During our nightly monitoring of M31 in the intermedi- 
ate Palomar Transient Factory (iPTF; hereafter called simply 
PTF; Law et al. 2009), we discovered a transient at the loca- 
tion of a known nova in M3 1 (Tang et al. 2013), and confirmed 
it to be a RN with a recurrence time of lyr with four optical 
outbursts detected in PTF from 2009 to 2013. During the op- 
tical novae, only a fraction of envelope is ejected (Starrfield 
et al. 1974), while the remaining envelope is expected to con- 
tinue hydrogen burning. As the ejected envelope expands, the 
ejecta becomes optically thin, and a supersoft source (SSS) 
powered by hydrogen burning is expected to emerge after the 
optical nova (Sala & Hernanz 2005; Wolf et al. 2013). Recent 
observational work has now made it clear that all novae have 
an extended supersoft phase whose duration and temperature 
depend solely on the WD Mass (Pietsch et al. 2005; Orio et al. 
2010; Henze et al. 2013c). We show here that this apparent 

14 http : //www . mpe . mpg . de/ { ~ { } }m31novae/opt/m31/ 

M31 table.html 


2 


Tang et al. 


transient behavior is an excellent match to the supersoft phase 
of a recurrent nova, and that the optical novae must have been 
missed in 1992, 1993 and 2001 when X-ray outbursts were 
seen (White et al. 1995; Williams et al. 2004). We describe 
the PTF discovery and optical follow-up observations in §2. 
Archival optical and X-ray studies are presented in §3 and §4, 
respectively. Swift observations and analyses are presented in 
§5. Theoretical modeling is presented in §6. Our conclusion 
is in §7. 

2. PTF DISCOVERY AND OPTICAL FOLLOW-UP 
OBSERVATIONS 

On 2013 Nov 27.08 UT 15 , we detected a transient at a = 
00''45"'28.89 J , 6 = 41°54'10"2 with R =18.9 mag (Tang et al. 
2013) in the nightly monitoring of M31 in the PTE using the 
48-inch telescope at Palomar. It brightened to R =18.3 mag 
on Nov 28.08. No source was detected at the same loca- 
tion to R <21 mag in PTF images taken on Nov 26.08 and 
Nov 25.29. There was no detection in 270 nightly PTF R- 
band images taken between 2013 May 19 to Nov 13 to a 
similar depth during non-bright time. The transient is co- 
incident within measurement uncertainties with the reported 
positions of three optical novae or novae candidates, i.e. 
He/N nova M31N2012-10a (Nishiyama & Kabashima 2012; 
Shafter et al. 2012), nova candidates M31N2011-10e (Ko- 
rotkiy & Elenin 2010; Barsukova et al. 2011) and M3 IN 
2008-12a (Nishiyama & Kabashima 2008). It is also coin- 
cident with the position of a ROSAT recurrent supersoft tran- 
sient RX J0045. 4+4154 (White et al. 1995), which was the 
first discovery of outbursts from this source. Therefore, we 
refer to the transient as RX J0045.4+4154 hereafter. 

Following our discovery, we initiated rapid photometric and 
spectroscopic follow-up observations. We obtained Bg'r'i' 
observations on the Palomar 60-inch telescope (P60; Cenko 
et al. 2006) on Nov 29, Nov 30, and Dec 1, and VI images on 
the Low Resolution Imaging Spectrometer (LRIS; Oke et al. 
1995) mounted on the Keck I 10 m telescope on Dec 4. Final 
reduction of the PTF images was performed using a forced- 
position Point-Spread Function (PSF) photometry pipeline 
(Ofek et al. 2012, Masci et al. in prep.). The P60 and LRIS 
images were reduced using aperture photometry. We cali- 
brated the LRIS images using the Local Group Survey catalog 
(LGS; Massey et al. 2006), and obtained a refined position for 
the transient of a = 00'‘45 m 28.847\ 6 = 41°54'10"05 with un- 
certainty of 0.1" (dominated by systematic uncertainties; the 
scattering of positions on the 3 LRIS images are « 0.01"). 
The optical light curve in 2013 is shown in the top panel in 
Figure 1. 

Optical spectroscopic follow-up of RX J0045. 4+4154 was 
undertaken with the Deep Imaging Multi-Object Spectro- 
graph (DEIMOS; Faber et al. 2003) mounted on the Keck II 
lOmtelescope on Nov 29, 2013 (1.3 days post peak), the Kast 
double spectrograph (Miller & Stone 1993) on the Shane 3 m 
telescope at Lick Observatory on Dec 1 (3.2 days post peak), 
and LRIS mounted on the Keck 1 10 m telescope on Dec 2 (4. 1 
days post peak). The spectral resolution is 4.2 A, 6 A (blue 
side) to 1 1 A (red side), and 7 A, for the DEIMOS, Kast, and 
LRIS spectrum, respectively. Spectra were reduced with stan- 
dard IRAF routines. 

The spectroscopic series of RX J0045.4+4154 is shown in 
Figure 2. It showed strong emission lines of Balmer series, 
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He I, He II, and N III. The observed optical lines place it into 
the “He/N" class of Williams (1992). The spectra are simi- 
lar to the HET spectrum taken 0.6 days post peak of the 2012 
outburst (M31N2012-10a; Shafter et al. 2012). Such similar- 
ity suggests that RN outbursts only depend on system param- 
eters like the WD mass and binary properties. Compared with 
the DEIMOS spectrum. He I lines weakened significantly in 
the LRIS spectrum, suggesting a decrease in the shell ioniza- 
tion during this period. The FWHM of Her emission line de- 
creased, from 2600 km s _1 in the DEIMOS spectrum, to 1900 
km s -1 in the LRIS spectrum, suggesting decreased ejecta ve- 
locity. Such velocity is at the low end of He/N nova (Shafter 
et al. 2011), and is at the low end of RNe (see e.g. Kato & 
Hachisu 2003; Yamanaka et al. 2010). 

The galactic extinction along the line of sight of RX 
J0045. 4+4154 is A R = 0.134 and A ? = 0.205 (Schlafly & 
Finkbeiner 2011), which is the lower limit of extinction of the 
object. To estimate the line-of-sight extinction contribution 
from M31 from the LRIS spectrum, we subtract the contin- 
uum, and measure the flux ratio of Hcr/H/3 to be 3.89. As- 
suming case B recombination (optically thick in Lycr Lines, 
optically thin in all other Hydrogen lines), the extinction is 
A r = 0.65 ±0.23 andA g = 1.0 ±0.35 (uncertainty comes from 
the range in expected ratios for Case B of 2.76-3.30; Oster- 
brock & Ferland 2006), which corresponds to Nu= 1 .5 ± 0.5 x 
10 21 cm -2 (Predehl & Schmitt 1995). If the nebula is not op- 
tically thin in Balmer lines or if the collisional excitation is 
non-negligible, a higher Ha/H/3 ratio is expected relative to 
the case B recombination, leading to an overestimated extinc- 
tion. Therefore, the above dust extinction and hydrogen col- 
umn density should be regarded as upper bounds. 

3. ARCHIVAL OPTICAL OBSERVATIONS 
3.1. Past Outbursts 

Archival work within PTF revealed other outbursts. We 
confirmed outbursts in Oct 2012 (He/N nova M31N2012- 
10a; Nishiyama & Kabashima 2012; Shafter et al. 2012) 
and Oct 2011 (nova candidate M31N201 l-10e; Korotkiy & 
Elenin 2010; Barsukova et al. 2011), and identified another 
un-reported outburst in Dec 2009. The transient is also coin- 
cident with the reported positions of nova candidate M3 IN 
2008- 12a (Nishiyama & Kabashima 2008). We also went 
back to Arp’s survey (1956) of M31 novae in the 1950s and 
found nothing at this location in his catalog. 

The observed outbursts are listed in Table 1. The recur- 
rence time is 341 days from ROSAT (White et al. 1995), and 
the recurrence times between the five most recent outbursts 
(2008-2013) are 342, 689, 363, and 404 days, respectively. 
Assuming a recurrence time between 330-410 days, we ex- 
pect an outburst in 2010 during Oct 29 to Nov 27. During 
this time, there are 30 PTF images (R band, typical limiting 
mag 20.5-21 mag), 6 images taken by P60 ( g ' and r' bands, 
typical limiting mag 21 -22 mag), and 6 images taken by 
MegaCam on the Canada-France-Hawaii Telescope (CFHT; 
Boulade et al. 2003) on Oct 31, 2010 in u' , g' , and r' with 
200s to 600s exposures (limiting mag ss 24 mag). No outburst 
was detected. If there is an outburst with light curve similar 
to the 2013 and 2011 novae, it will be detected by PTF during 
t—to = -1 to 1 (where to is the optical peak time), by P60 dur- 
ing t —to = — 1 to 1, and by CFHT/MegaCam during 1 — to = -1 
to 7. Only two possible windows are left if there was a nova 
in 2010, i.e. Nov 10-12, or Nov 20-27. 

The optical light curves of the 2009, 2011, 2012, and 2013 
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outbursts are shown in Figure 1, and are similar. We have 
daily sampled g', and R band light curves covering the rise, 
peak, and decay of the 2011 and 2013 outbursts. The de- 
rived ?2 (time to decay from the peak by 2 mag) is t 2 , g =2.0 
days in the 2011 outburst, and f 2 j? =2.1 days in the 2013 out- 
burst, making it among the fastest novae (Cao et al. 2012). 
The declines within 2 mag from the peak are more or less lin- 
ear, and thus linear regression is used to derive to. The peak 
magnitude in the 2013 outburst is R = 18.34 ± 0.08 mag at 
MJD = 56624.076 from PTF data. The peak magnitudes in the 
2011 outburst is g = 18.51 ±0.09 mag at MJD = 55857.137 
from PTF data, and R = 18.18 ±0.08 mag at MJD = 55857.121 
from Barsukova et al. (2011). 

The simultaneous or nearly simultaneous photometric mea- 
surements in the 2011 and 2013 outbursts are listed in Table 
1. Assuming an intrinsic effective temperature of T k, 8200 K 
as observed in typical novae, the inferred extinction is A R « 
0-0.7 (V H < 1.6 x 10 21 cm -2 ). 

Adopting an extinction of A R = 0.45 and A g = 0.7 (/V H = 
1.0 x 10 21 cm -2 ), which is consistent with both the Balmer 
decrement and available colors of the outbursts, and given the 
distance modulus of M31 of 24.4 (Vilardell et al. 2010), the 
peak absolute magnitude of the 2013 and 2011 outbursts is 
then M r = -6.5 mag, and M g = -6.6 mag, respectively. With 
the fastest decline rate (1 mag d -1 ) yet faint peak magnitude 
among novae, RX J0045. 4+4154 is an extreme outlier in the 
canonical maximum-magnitude-rate-of-decline (MMRD) re- 
lation (della Valle & Livio 1995). As discussed by Kasliwal 
et al. (2011), such faint and fast novae can arise from pro- 
genitors containing high accretion rate and relatively massive 
white dwarfs (thus lower envelope masses at Hydrogen igni- 
tion; see e.g. Wolf et al. 2013), as we expect for recurrent 
novae. 

3.2. Archival HST Observations 

Following the candidate progenitor reported in Williams 
et al. (2013), we re-analysed the Hubble Space Telescope 
(HST) archival exposures to investigate photometric prop- 
erties of RX J0045. 4+4154 in quiescence. Besides the 
ACS/WFC data (Filters: F475W and F814W; Date; 2010 
August 7) mentioned in the report, we made use of newer 
ACS/WFC (Filters; F475W and F814W; Date: 2012 Jan- 
uary 10) and WFC3/UVIS&IR (Filters: F110W, F160W, 
F275W, and F336W; Date: 2011 January 25 and August 
31) observations in the analysis. To search for the progeni- 
tor, accurate astrometry with precision down to 0.02"-0.07" 
between our Keck/LRIS image (see §2) and the HST im- 
ages were performed based on 5-18 bright reference stars. 
Within the Keck/LRIS error circle (astrometric uncertainty 
dominated) of RX J0045.4+4154, a clear source is present 
in the optical/UV-band HST images but is undetectable in 
the IR band. We measured magnitudes of the detected 
source through the Dolphot PSF photometry package (Dol- 
phin 2000) with HST-dedicated parameters suggested in the 
manual. All calculated photometric measures are listed 
in Table 3. The mean magnitudes for the optical/UV fil- 
ters are F275W=23.07, F336W=22.96, F475W=24.24, and 
F814W=23.91 16 . For the IR bands, we estimated the upper 
limits by examining the faintest stars detected by the HST ex- 
posures, which gives us F110W > 24.44 (exposure: 800s) 
and F160W > 25.22 (exposure: 1700s). Variability is present 

16 See Table 3 for details. 


in all optical/UV bands with amplitudes up to 0.3 mag for 
F275W, 0.6 mag for F336W, 0.5 mag for F475W, and 0.2 mag 
for F814W, likely from the high accretion rate disk. There 
is a significant change in color between 2010 and 2012 with 
A(B-I) « -0.3 suggesting the system was possibly in a dif- 
ferent phase of the nova recurrent cycle or/and at a different 
accretion rate. The results are consistent with an independent 
study on the same object by Darnley et al. (2014). 

4. ARCHIVAL X-RAY DATA 

In the M31 second ROSAT PSPC survey catalog (Sup- 
per et al. 2001), a recurrent supersoft X-ray transient RXJ 
0045.4+4154, identified by White et al. (1995), is spatially 
coincident with the optical RN (offset: 4"; uncertainty 13"). 
RXJ 0045.4+4154 showed two distinct supersoft X-ray out- 
bursts of duration « 5 days in February 1992 and January 
1993 (340 days apart) (see White et al. 1995, figure 2), 
in which rapid luminosity changes over the 5-day intervals 
were clearly seen while a comparably bright X-ray source 
was completely absent before or after the outbursts. White 
et al. (1995) interpreted the supersoft X-rays as the conse- 
quence of thermonuclear burning of hydrogen on the surface 
of a WD and estimated the WD atmosphere temperature to 
be k7bb ~ 90 eV by fitting the spectrum with a blackbody 
model, inferring a peak luminosity (absorption corrected) 
of L x ~ 10 38 erg s -1 (0.1-2 keV) and a blackbody radius of 
Rbb ~ 4.7 x 10 8 cm (/J Mtl = 780 kpc assumed). From an O 
VIII ionization edge detection at 0.87 keV, RXJ 0045.4+4154 
was believed to be a massive WD with M > 1.3 Mq (i.e., 
log(g) > 8.75). 

Williams et al. (2004) reported a Chandra ! HRC transient 
detection in Sep 2001, in which the X-ray transient, named 
nl-85, was mostly quiescent ( L x < 2 x 10 36 erg s -1 ), but was 
detected at one epoch with L x « 6 x 10 37 erg s -1 (see Williams 
et al. 2004, Figure 6). Their reported coordinates of nl-85 
is 9" offset from the optical position of RX J0045. 4+4154, 
seemingly making an association with RX J0045. 4+4154 un- 
likely. However, we registered the Chandra ! HRC image to 
the Swift/XRT WCS frame to investigate possible association. 
Using four common X-ray sources detected by both the in- 
struments in the field, the astrometric corrected offset of nl-85 
and RXJ 0045.4+4154 dropped to about 1", which is smaller 
than the positional uncertainty of the XRT position (1.6"), 
hence, we are confident of the association. 

5. SWIFT OBSERVATIONS 
5.1. Swift XRT and UVOT light curves 

Given the short X-ray outburst duration observed in ROSAT 
(White et al. 1995), we launched high cadence ToO X-ray 
and UV observations with the Swift observatory (Gehrels et al. 

2004) . Another Swift campaign was carried out by M. Henze 
and collaborators (Henze et al. 2013a) following our optical 
discovery (Tang et al. 2013). RX J0045. 4+4154 was observed 
with the Swift X-ray telescope (XRT; Burrows et al. 2005) 
and the Ultraviolet/Optical Telescope (UVOT; Roming et al. 

2005) in a series of 36 observations beginning on 2013 Dec 3 
(5 days after the optical peak). Typically 4-6 ks observations 
were taken in 1-2 ks snapshots each day from Dec 3 to Dec 
16, and 7-8 ks observations once every 3 days from Dec 17 
to Dec 23. The total exposure time is 97.6 ks. 

A highly variable source is detected in both XRT and UVOT 
at the optical position of RX J0045. 4+4154 (see also Henze 
et al. 2013b). The enhanced position given by the online 
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Swift XRT pipeline software of the UK Swift Science Data 
Centre at the University of Leicester (Evans et al. 2007) is 
a = 00 / ’45 m 28.80\ S = 41°54'08"8 with a 90% error of 1.6", 
which is consistent with the optical position of the nova. The 
Swift XRT light curve of RX J0045. 4+4154 given by the same 
pipeline is shown in black in the top panel of Figure 3. The 
first detection in optical is 1 day before the optical peak, and 
we have a 5<r upper limit of R > 21 mag 2 days before the 
optical peak. Hence the onset of the thermonuclear runaway 
(TNR) is 1-2 days before the optical peak. Therefore, from 
the XRT light curve, we measured the turn-on time (defined 
as the time it takes for the source to emerge as a SSS after 
the onset of TNR) of <6-7 days, and turn-off time (defined 
as the time it takes for the source to disappear in X-rays after 
the onset of TNR) of 18-19 days. Both timescales are among 
the shortest ones measured for novae (Henze et al. 2013c). 

The UVOT data were reduced with the HEASoft V6.15 
package and with the calibration files released in January 
2013. Aperture photometry was performed using UVOT- 
PRODUCT, with a 5" radius for the source, and an annulus 
with inner radius 27" and outer radius 35" for background, as 
recommended by Poole et al. (2008). No bright UV source is 
located in the background annulus. The resulting UVOT light 
curves are shown in the bottom panel of Figure 3. The magni- 
tudes are in the Vega system. It showed significant short-time 
variations, notably 1 mag variations on hourly timescales in 
the uvw2 filter (1928 ±657 A; Poole et al. 2008) during the 
5 snapshots in the first Swift observation. It also showed a 
relatively monotonical decline on longer timescales (a few 
days) after the first Swift observation. The long-term varia- 
tion amplitude during the Swift campaign is about 2 mag, 1 .4 
mag and 1 mag in the uvw2, uvm2 (2246 ±498 A) and uvwl 
(2600 ±693 A) filters, respectively. No significant variation is 
detected in the u band (3465 ± 785 A). The relatively smaller 
variation amplitudes in UV compared with X-ray are consis- 
tent with the scenario of decreasing effective temperature in 
the late stage of SSS. 

5.2. Effective temperature and column density evolution 

We obtained the Swift data including the corresponding 
ancillary files from the Swift quick look data archive, ex- 
tracted the level-2 event files to spectra using HEASoft 
version 6.14, and performed spectral fittings using XSPEC 
version 12.8.1 with the an absorbed blackbody model (i.e., 
phabs*bbobyrad). Given that the X-ray variability is high 
over the entire supersoft phase, we split the spectrum into 
parts according to the observing time to investigate the lu- 
minosity and temperature evolutions of the system. We di- 
vided the observations into three groups (Table 4) and fit- 
ted each with an absorbed blackbody. Due to a low photon 
count, we assumed two possible V H , he- 1.0 x 10 21 cm -2 and 
1 .5 x 10 21 cm -2 , to increase the degrees of freedom of the fits. 
The assumed /V H are consistent with the extinction measured 
from the Balmer decrement (N H < 1.5 x 10 21 cm -2 ) and op- 
tical colors ( Nh= 0-1.6 x 10 21 cm -2 ), as well as the results 
from an independent study on the same object by Henze et al. 
(2014) who estimated a Nh of 1.1 — 1.6 x 10 21 cm -2 from the 
Swift data. 

As listed in Table 4, the fits show significant variabilities 
in temperature and luminosity and all best-fit parameters are 
consistent with the ROSAT outbursts. The limited number of 
XRT photons and our lack of knowledge of the metallicity at 
the stellar photosphere inhibited our use of the more accurate 


stellar atmosphere models of Rauch et al. (2010) that were 
constructed for hot and massive WDs. 

6. HIGH M AND M NOVA MODELS 

The short nova recurrence time, rapid evolution as an X- 
ray source and high surface temperature during the SSS phase 
all consistently point to a high mass WD. To quantify just 
how large the WD mass must be to explain the observations, 
we undertook an expansion of the recent work of Wolf et al. 
(2013) using the Modules for Experiments in Astrophysics 
(MESA rev. 5596; Paxton et al. 2011, 2013). We simulated 
WDs with M = 1 .30, 1 .32, 1 .34, and 1 .36 M 0 accreting mate- 
rial with solar composition, focusing on models that yielded 
the observed recurrence time of 1 yr, yielding an accretion 
rate range of 1.7 x 10 -7 < M/M 0 yr -1 < 3.3 x 10 -7 . These 
model WDs have core temperatures of T c = 3 x 10 7 K, ex- 
cept for the 1.36 M 0 model, which had T c = 6 x 10 7 K. The 
value of T c does not impact the outcome at these high M’s 
due to the heat buffer created by the even hotter helium layer 
(r He ~ 10 s K) (Prialnik & Kovetz 1995; Cassisi et al. 1998; 
Piersanti et al. 1999; Yaron et al. 2005; Wolf et al. 2013). The 
mass loss prescription during the novae events is a super Ed- 
dington wind, as described in Denissenkov et al. (2013) and 
Wolf et al. (2013). 

These calculations immediately show that the WD in RX 
J0045. 4+4154 has a mass of at least 1.3 M 0 , as lower mass 
WDs could not yield a minimum recurrence time as short as 
1 yr. The M’s for the 1.32, 1.34, and 1.36 M 0 models that 
yielded f recur = 1 yr were, respectively, 3.1 x 10 -7 , 2.1 x 10 -7 , 
and 1.7 x 10 -7 M 0 yr -1 . These models retained ss 30% of 
the accreted material through the outburst, yielding effective 
accretion rates onto the helium layer of Mwd ~9x 10 -8 , 9 x 
10 -8 , and 6 x 10 -8 M 0 yr -1 . Assuming this matter ultimately 
stays on the WD and that the effective accretion rate remains 
constant, these models would evolve to M = 1.37 M 0 within 
5 x 10 5 years. 

We next studied the location of these high mass WDs in 
the HR diagram during the SSS phase, all of which enter 
the SSS phase after only about 10-20 days from the onset 
of the TNR. As listed in Table 4, the measured 7 dT from 
the two absorbed blackbody models are consistent with each 
other within uncertainties, though there is a modest discrep- 
ancy in the blackbody luminosities. As shown in Figure 5, 
all theoretical tracks are consistent with the measured 7 0 f 
of RX J0045.4+4154, with either N H = L0 x 10 21 cm -2 or 
A'n= 1.5 x 10 21 cm -2 . However, the WDs have different rates 
of evolution in the HR diagram depending on their mass, as 
denoted by the solid points that denote the elapsed time be- 
tween each point. For example, the 1.34 M 0 WD (orange 
dotted line and triangles) evolves through the SSS phase in 
about 15 days, whereas the 1.36 M 0 WD is in a SSS phase 
for only 8 days. Clearly, the 1.30 M 0 WD evolves far too 
slowly. So, this comparison again points to WD masses of 
1.30 < M/Mq < 1.36 in our super Eddington wind model, 
with mass loss due to Roche lobe overflow indicating more 
massive WDs (Wolf et al. 2013). The turnoff times could be 
shorter if there was substantial mixing of elements heavier 
than hydrogen (e.g. He, C, O) into the burning layer during 
the TNR (Sala & Hernanz 2005), pointing to WD masses at 
the lower end, i.e. M = 1 .30 M 0 . 

7. CONCLUSION 
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As a recurrent nova, RX J0045. 4+4154 has many remark- 
able features, with the shortest recurrence time of 1 year, rapid 
turn-on and turn-off of the stable burning SSS phase, and the 
highest peak effective temperature during the SSS (100-110 
eV). We showed here that these are all remarkably consistent 
with theoretical models of hydrogen thermonuclear flashes on 
a WD with a mass in the range of 1.30 < M/M@ < 1.36. Se- 
curely identifying such a massive WD is key to our under- 
standing of the larger problem of deciding how some massive 
WDs lead to explosions as Type la SNe (if the core is car- 
bon rich) or undergo an accretion-induced collapse (when the 
core is composed of O/Ne). At the accretion rate inferred 
from our theoretical calculations, and assuming (from the du- 
ration of the supersoft phase) that about 30% of the accreted 
material stays on the WD, it will take less than a million years 
for the core density of the accreting WD to reach values ade- 
quate for either an unstable carbon ignition or onset of elec- 
tron capture (if an O/Ne core). This assumes that mass loss in 
the inevitable intervening unstable helium flashes is negligible 
(Kato & Hachisu 2004). 

This remarkable system will certainly undergo additional 
outbursts, and we hope that the next one (likely in Nov-Dec 
2014) will be studied in even more detail, especially in the X- 
rays where spectra can certainly reveal much more about the 
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WD mass and surface composition. More detailed theoreti- 
cal modeling is certainly justified, including simulations that 
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This work was supported by the National Science Founda- 
tion under grants PHY 11-25915, AST 11-09174, and AST 
12-05574. Most of the MESA simulations for this work were 
made possible by the Triton Resource. The Triton Resource 
is a high-performance research computing system operated by 
the San Diego Supercomputer Center at UC San Diego. This 
research used resources of the National Energy Research Sci- 
entific Computing Center, which is supported by the Office of 
Science of the U.S. Department of Energy under Contract No. 
DE-AC02-05CH1 1231. AKHK is supported by the National 
Science Council of the Republic of China (Taiwan) through 
grant NSC101-21 19-M-008-007-MY3. MMK acknowledges 
generous support from the Hubble Fellowship and Carnegie- 
Princeton Fellowship. We are grateful to the Swift Team for 
the superb timely scheduling of the observations and provid- 
ing data and analysis tools, and to Bill Paxton for his devel- 
opment of MESA. 

Facilities: PO:1.2m (PTF), 1.5m, Keck (DEIMOS, LRIS), 
Swift (XRT, UVOT). 


Nishiyama, K., & Kabashima, F. 2008, Central Bureau for Astronomical 
Telegrams, http : / / www . cbat . eps . harvard . edu/ iau/CBAT. 
M31 ,html#2008-12a 

— . 2012, Central Bureau for Astronomical Telegrams, 

http : //www . cbat . eps . harvard. edu/unconf/ followups/ 
JO 04 52 88 4+4 15 4 095. html 
Ofek, E. O., Laher, R„ Law, N„ et al. 2012, PASP, 124, 62 

Oke, J. B„ Cohen, J. G„ Carr, M„ et al. 1995, PASP, 107, 375 
Orio, M„ Nelson, T., Bianchini, A., Di Mille, F., & Harbeck, D. 2010, ApJ, 
717,739 

Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of gaseous nebulae 
and active galactic nuclei 

Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 
Paxton, B., Cantiello, M., Arras, P, et al. 2013, ApJS, 208, 4 
Piersanti, L., Cassisi, S., Iben, I. J.. & Tornambe, A. 1999, ApJ, 521, L59 
Pietsch, W„ Fliri, J., Freyberg, M. J., et al. 2005, A&A, 442, 879 
Poole, T. S„ Breeveld, A. A., Page, M. J., et al. 2008, MNRAS, 383, 627 
Predehl, P„ & Schmitt, J. H. M. M. 1995, A&A, 293, 889 
Prialnik, D„ & Kovetz, A. 1995, ApJ, 445, 789 

Rauch, T., Ringat, E., & Werner, K. 2010, ArXiv e-prints, arXiv:1011.3628 
Roming, P. W. A., Kennedy, T. E., Mason, K. O., et al. 2005, 

Space Sci. Rev., 120, 95 
Sala, G„ & Hemanz, M. 2005, A&A, 439, 1061 
Schaefer, B. E. 2010, ApJS, 187, 275 
Schlafly, E. F„ & Finkbeiner, D. P. 201 1, ApJ, 737, 103 
Shafter, A. W., Curtin, C., Pritchet, C. J., Bode, M. F., & Darnley, M. J. 

2013, ArXiv e-prints, arXiv:1307.2296 
Shafter, A. W., Hornoch, K., Ciardullo, J. V. R.. Darnley, M. J., & Bode, 

M. F. 2012, The Astronomer’s Telegram, 4503, 1 
Shafter, A. W., Darnley, M. J., Hornoch, K., et al. 2011, ApJ, 734, 12 
Starrfield, S„ Sparks, W. M„ & Shaviv, G. 1988, ApJ, 325, L35 
Starrfield, S„ Sparks, W. M„ & Truran, J. W. 1974, ApJS, 28, 247 
Supper, R., Hasinger, G., Lewin, W. H. G., et al. 2001, A&A, 373, 63 
Tang, S., Cao, Y., & Kasliwal, M. M. 2013. The Astronomer’s Telegram, 
5607, 1 

Vilardell, F., Ribas, I., Jordi, C., Fitzpatrick, E. L., & Guinan, E. F. 2010, 
A&A, 509, A70 

White, N. E., Giommi. P., Heise, J., Angelini, L., & Fantasia, S. 1995, ApJ. 
445 LI 25 

Williams, B. F„ Garcia, M. R„ Kong, A. K. H„ et al. 2004, ApJ, 609, 735 
Williams, R. E. 1992, AJ, 104, 725 

Williams, S. C., Darnley, M. J., & Bode, M. F. and Shafter, A. W. 2013, The 
Astronomer’s Telegram, 5611, 1 

Wolf, W. M„ Bildsten, L.. Brooks, J., & Paxton, B. 2013, ApJ, 777, 136 
Yamanaka, M., Uemura, M., Kawabata, K. S., et al. 2010, PASJ, 62, L37 
Yaron, O., Prialnik, D„ Shara, M. M., & Kovetz, A. 2005, ApJ, 623, 398 


6 


Tang et al. 





s 

s s 


18 


9 

1 

1 '1 

2013 outburst 

19 


• 

e ! . 
+ 


f 0 =56624.076 

Di 

ro 20 



i 

t 2R = 2.13 d 

21 


l 



• 

9 

22 






18 


* * 



2012-10a 

19 




t 0 =56219.730 - 






ro 20 





T 





s 

* 1 T T T I 

21 





1 J 1 ^ 

22 






18 


•s 



2011-10e " 

™ 19 
CD 

fD 20 


% 


t 0 = 55857.137 - 

t 2 , g = 1-97 d . 




S’ 



21 

" T T 



% 

* I l ^ 

22 

. 1 i 




i 1 

18 





2009 outburst 

19 


• 



f 0 =55168.174 - 

cn 





ro 20 





- 



T 




21 


T 



- 






i 

22 





- 


-4 -2 0 2 4 6 8 10 

Days after optical peak 


PTF 

P60 

LRIS 

Others 

9 

r 

i 

B 

V 

R 

I 

clear 


Fig. 1. — Optical light curve of RX J0045.4+4154. Four outbursts were detected in PTF from 2009 to 2013. PTF data are shown in circles (green for g ' band, 
and red for R band). P60 follow-ups are shown in diamonds (blue for B band, green for g' band, red for r' band, and magenta for i' band). LRIS follow-ups are 
shown in hexagons (cyan for V band, and magenta for I band). Published photometry by others (Korotkiy & Elenin 2010; Barsukova et al. 2011; Nishiyama & 
Kabashima 2012; Shafter et al. 2012; Williams et al. 2013) are shown in squares (black for clear, blue for B band, cyan for V band, red for R band, and magenta 
for I/i' band). Red arrows are 5cr uplimits from PTF R band images. On the top axis, the epochs of spectroscopic follow-up of the 2013 outburst are indicated 
by “S”. 


TABLE 1 

List of observed outbursts from RX J0045. 4+4154. 


t(),opt a 

(UT) 

fo^ b 

(UT) 

Time since last observed outburst Source 

(days) 

Reference 


2013-12-05 


X-ray (Swift) 

this paper; Henze et al. (2014) 

2013-11-28 


404 

Optical (PTF) 

Tang et al. (2013); this paper 

2012-10-19 


363 

Optical 

Nishiyama & Kabashima (2012); Shafter et al. (2012) 

2011-10-23 


689 

Optical 

Korotkiy & Elenin (2010); Barsukova et al. (201 1); this paper 

2009-12-03 


342 

Optical (PTF) 

this paper 

2008-12-26 



Optical 

Nishiyama & Kabashima (2008) 


2001-09-08 


X-ray ( Chandra ) 

Williams et al. (2004) 


1993-01-11 

341 

X-ray ( ROSAT) 

White etal. (1995) 


1992-02-05 


X-ray ( ROSAT) 

White etal. (1995) 


“Time of the optical peak. 
b Time of the X-ray peak. 
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Fig. 2. — Optical spectra of RX J0045.4+4154 during the 2013 outburst. Left panels: arbitrarily normalized spectra. Each spectrum is labelled with the 
observation date, the telescope and the instrument. Right panels: Ha and H (3 line profiles. The Kast spectrum is offseted by +0.5, and the DEIMOS spectrum is 
offseted by +1.0. The zero levels of the flux are shown in horizontal dotted lines. The position of zero velocity is marked by the vertical dashed line. 


TABLE 2 

Optical photometry of RX J0045.4+4154 during the 2013 and 201 1 outbursts with (nearly) simultaneous multi-band observations. 


MJD 

t-to 

(days) 

B 

(mag) 

V 

(mag) 

R I 

(mag) (mag) 

g’ 

(mag) 

r’ 

(mag) 

i’ 

(mag) 

Ref. “ 

2013, t 0 = 

56624.08 in MJD 







56624.93 

0.85 

19.61 ±0.01 

19.65 ±0.02 




19.29 ±0.02 

W13 

56625.26 

1.18 

20. 15 ±0.20 



19.82± 0.18 

19.54±0.16 

19.47 ± 0.17 

P60 

56626.15 

2.07 

20.67 ±0.10 


20.32± 0.16 

20.67 ± 0.07 



PTF and P60 

56627.32 

3.24 

20.63 ±0.17 



20.62 ±0.10 



P60 

56630.34 

6.26 


2 1.30 ±0.02 

20.92 ±0.05 




LRIS 

2011, to = 

55857.1 in MJD 







55857.1 

0 



18.18 ±0.08 

18.51 ±0.09 



Bll; PTF 

“References: W13 

is Williams et al. 

2013; Bll is 

Barsukova et al. 2011. Others 

are either PTF photometry or 

our follow-up using Palomar 60- 


Keckl/LRIS. 
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Time after optical peak (days) 


Fig. 3. — Top: Swift XRT light curve of RX J0045.4+4154. Arrows are 3cr upper limits. Vertical dotted lines mark the time bins used in spectral analysis 
(Figure 4 and Table 4). Bottom: Swift UVOT light curves of RX J0045.4+4154. The uvm2, uvwl, and u band light curves are shifted by constant values as 
indicated in the legend. Open symbols with arrows are 3cr upper limits. 
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Fig. 4. — The SwiftfXKT spectra (from 0.3 to around 1 keV) of RX J0045. 4+4154 taken at different time bins during the supersoft phase are shown with their 
best-fit black body spectral models assuming a fixed Nu=l.5 x 10 21 cm -2 . Black for days 5-8 after the optical peak, red for days 8-11, and green for days 1 1-17. 
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Fig. 5. — Comparison of observations vs theoretical models on the Luminosity (bolometric) vs ksT e ff plane during the SS phase. Best-fit results from the Swift 
XRT data using an absorbed blackbody model with fixed Nr=1.0 x 10 21 cm -2 and Nu=l.5 x 10 21 cm -2 are shown in black upside down triangles and red crosses, 
respectively, with the top-right point from the early stage (5-8 days after optical peak), the middle-left point from the middle stage (8-1 1 days after optical peak), 
and the lower-right point from the late stage (11-17 days after optical peak). Note that Lx(03— 1 keV) in Table 4 are converted to bolometric luminosity assuming 
blackbody radiation with best-fit temperatures. Evolutionary tracks for the four MESA models with r re cur ~ 1 yr are shown in thick colored lines with markers; 
all the tracks evlove from higher luminosity to lower luminosity. Different lines indicate different WD mass and M combinations (indicated in the legend in units 
of Mq and Mq yr -1 ), and the markers separate periods of equal time which are also indicated in the legend. More massive stars yield shorter SSS phases due to 
the smaller envelope masses as well as the more vigorous burning (Wolf et al. 2013). Also shown in gray are lines of constant radius. 
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TABLE 3 

HST Observations of RX J0045. 4+4154 


Filter ID 

Observing Time 

Days after“ 

Days before 3 

Exposure 

Vega Magnitude 


(yyyy-mm-dd hh:mm:ss) 

previous nova 

next nova 

(s) 

(mag) 

F275W 

2011-01-25 04:56:24 

418 (74) 

271 

350 

23.244± 0.1 18 


2011-01-25 05:20:52 

418 (74) 

271 

660 

23.078 ±0.072 


2011-08-31 12:18:37 

636 (292) 

53 

350 

23.075 ±0.103 


2011-08-31 12:41:25 

636 (292) 

53 

575 

22.955 ±0.079 

F336W 

2011-01-25 04:44:43 

418 (74) 

271 

550 

23.309 ±0.054 


2011-01-25 05:04:54 

418 (74) 

271 

800 

23. 118 ±0.040 


2011-08-31 12:06:56 

636 (292) 

53 

550 

22.715 ±0.042 


2011-08-31 12:27:07 

636 (292) 

53 

700 

22.738 ±0.037 

F475W 

2010-08-07 12:27:39 

247 

431 (87) 

600 

24.041 ±0.022 


2010-08-07 12:40:20 

247 

431 (87) 

370 

24.041 ±0.029 


2010-08-07 12:49:08 

247 

431 (87) 

370 

24.035 ±0.029 


2010-08-07 12:57:56 

247 

431 (87) 

370 

23.989 ±0.028 


2012-01-10 02:45:00 

79 

284 

700 

24.448 ±0.028 


2012-01-10 02:59:21 

79 

284 

360 

24.451 ±0.038 


2012-01-10 03:07:59 

79 

284 

360 

24.410 ±0.037 


2012-01-10 03:16:37 

79 

284 

470 

24.478 ±0.033 

F814W 

2010-08-07 10:44:08 

247 

431 

350 

23.846 ±0.045 


2010-08-07 10:52:38 

247 

431 (87) 

700 

23.871 ±0.031 


2010-08-07 11:06:56 

247 

431 (87) 

455 

23.815±0.037 


2012-01-10 00:23:51 

79 

284 

350 

23.981 ±0.049 


2012-01-10 01:09:03 

79 

284 

800 

23.967 ±0.039 


2012-01-10 01:25:01 

79 

284 

550 

23.964 ±0.038 


“Days in parenthesis are the ones assuming there was a missing nova on Nov 12, 2010, which is the middle point between the 2009 and the 2011 novae. 


TABLE 4 

X-RAY SPECTRAL FITTING OF RX J0045. 4+4154 

Time a Nh Temperature (7bb) b Radius c {Rbb) Blackbody emission 11 (Lx) xt (dof) 
(days) (10 21 cnT 2 ) (eV) (10 8 cm) (10 38 ergs -1 ) (10 38 ergs _1 ) 


Single Absorbed Blackbody with a fixed Ah=l .0 x 10 21 cm 


5-8 

943 

4.7 3? 

i n+0.1 

a - z -0.3 


8_11 1.0 (fixed) 

112 + J, 

2-9 1[ 9 6 

1-23:2 

1.40 (42) 

11-17 

9571 

3 . 43:9 

o.73! 



Single Absorbed Blackbody with a fixed Ah=l .5 x 10 21 cm 


5-8 

00 
<1 

1 + 

7 T * -5 - 4 
' • 3 - 4.3 

? f) +0 ' 2 

Z.U_o.4 


8-11 1.5 (fixed) 

104^ 

4 ^ +2 - 9 
4 - 3 - 2.4 

1 *®-0.3 

1.27 (42) 

11-17 

OO 

00 

1 + 

C 1+3.5 

D. 1-2.9 

1 1+ 01 
1J 42 



“Time after the optical peak in 2013. 

b All uncertainty ranges are 90% confidence interval. 

c The radii of the blackbody were calculated using the best-lit normalizations (i.e., oc Rbb/^M 3 l Mth AvDl = 280 kpc). 

d The luminosities (0.3-1 keV) are absorption corrected with an assumption of Z>M 3 i = 780 kpc while the uncertainties were estimated by approximating 
ALunabs/T una bs — - ^Lbs/Lbs- 


